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*— i Abstract 

Systems biology uses large networks of biochemical reactions to model the functioning of biological 
£ — cells from the molecular to the cellular scale. The dynamics of dissipative reaction networks with many 

well separated time scales can be described as a sequence of successive equilibrations of different subsets of 
variables of the system. Polynomial systems with separation are equilibrated when at least two monomials, 
of opposite signs, have the same order of magnitude and dominate the others. These equilibrations and 
the corresponding truncated dynamics, obtained by eliminating the dominated terms, find a natural 
formulation in tropical analysis and can be used for model reduction. 
q Keywords: Tropical analysis, asymptotic analysis, chemical kinetics, systems biology. 

AMS subjects: Primary 14T05, 92C40, 92C42; Secondary 14M25. 

i — i 1 Introduction. 

Systems biology develops biochemical dynamic models of various cellular processes such as signalling, 
^{■^ metabolism, gene regulation. These models can reproduce complex spatial and temporal dynamic behavior 

observed in molecular biology experiments. The dynamics of multiscale, dissipative, large biochemical 
Q\ models, can be reduced to that of simpler models, that were called dominant subsystems [RGZL08, 

IGRZ101 IGR08| . Simplified, dominant subsystems contain less parameters and are more easy to analyze. 
£f) The notion of dominance is asymptotic and a natural mathematical framework to capture multiple 

asymptotic relations is the tropical analysis. Motivated by applications in mathematical physics [LM96 , 
systems of polynomial equations [StuQ2] . etc., tropical analysis uses a change of scale to transform nonlinear 
systems into piecewise linear systems. 

In this paper we provide some mathematical justifications for possible applications of the idea of 
tropicalization to systems biology models. 



2 Tropicalized chemical kinetics 

In chemical kinetics, the reagent concentrations satisfy ordinary differential equations: 

dx ■ 

= Fi(x), 1 < * < n. (1) 
Rather generally, the rates are rational functions of the concentrations and read 

Fi(x) = Pi(x)/Qi(x), (2) 
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where Pi(x) = ^2 aeA . di. a x a , Qi(x) = X)/3es bi,pxP , are multivariate polynomials. Here x a = x^x^ 2 ■ ■ ■ x" n , 

x^ = x^ 1 x^ 2 . . . x^ n , ai. a , bi^p, are nonzero real numbers, and Ai, Bi are finite subsets of W l . 
The special case of mass action kinetics is represented by 

F l {x)=P+{x)-Pr(x), (3) 

where P^(x), Pf{x) are positive coefficients polynomials, P i (x) = J2 a eA ± a fa xa > a ia > 0> and A t 
are finite subsets of N™. 

In multiscale biochemical systems, the various monomials defining reaction rates have different orders, 
and at a given time, there is only one or a few dominating terms. Therefore, it could make sense to replace 
polynomials with positive real coefficients X^qgA a aX a , by max-plus polynomials exp(max Q , 6y i(Zo(7(a a )+ < 
log(x),a >)). 

This heuristic can be used to associate a piecewise-smooth model to the system of rational ODEs ([T]), 
in two different ways. 

The first method was proposed in [NGVR12J and can be applied to any rational ODE system defined 

by 0,©: 

Definition 2.1. We call complete tropicalization of the smooth ODE system (JTj) , ([2]) the following piecewise- 
smooth system: 

dx ■ 

= DomPi(x)/DomQi{x), (4) 

where Dom (J^aeA a i,aX a ) = sign(ai tQ . max ) exp(m&x a( zA i (log(\a.i, a \)+ < u,a >)). Here u = (logxi, . . . , logx n ), 
< u,a > denotes the dot product, and a>i,a max > oc max € Ai denotes the coefficient of the monomial for 
which the maximum is attained. In simple words, Dom renders the monomial of largest absolute value, 
with its sign. 

The second method, proposed in jSCF + 09) . applies to the systems ([!]), ([3|. 

Definition 2.2. We call two terms tropicalization of the smooth ODE system ([l]),(|3| the following 
piecewise-smooth system: 

= DomP+(x) - DomPr(x), (5) 

The two-terms tropicalization was used in |SCF + 09] to analyse the dependence of steady states on the 
model parameters. The complete tropicalization was used for the study of the model dynamics and for 
the model reduction |N0VM2j |MMl2j. 

For both tropicalization methods, for each occurrence of the Dom operator, one can introduce a tropical 
manifold, defined as the subset of K™ where the maximum in Dom is attained by at least two terms. For 
instance, for n — 2, such tropical manifold is made of points, segments connecting these points, and half- 
lines. The tropical manifolds in such an arrangement decompose the space into sectors, inside which one 
monomial dominates all the others in the definition of the reagent rates. The study of this arrangement 
give hints on the possible steady states and attractors, as well as on their bifurcations [NGVR12] , 



3 Tropical equilibration and permanence 

In the general case, the tropicalization heuristic is difficult to justify by rigorous estimates. However, this 
is possible when the polynomials defining the rhs of the ODE system have dominant monomials, much 
larger than the other monomials. 

To simplify, let us consider the class of polynomial systems, corresponding to mass action law chemical 
kinetics : 

f = X>*~> <« 



2 



where otij are multi-indices, Pij are rate constants. 

In order to introduce orders, we consider that coefficients are integer powers of a small positive 
parameter e: 

P ij (e) = e^P ij . (7) 

where P^ do not depend or are 0(1) on e. 

We also suppose that the cone R> = {x : Xj > 0} is invariant under dynamics ([6| and initial data are 
positive: 

2,(0) > 6 > 0. 

The terms Q can have different signs, the ones with P^ > are production terms, and those with P^ < 
are degradation terms. 

From the biochemical point of view, the choice ^ is justified by the fact that biochemical processes 
have many, well separated concentration and time scales. The orders of different monomials defining the 
system ^ are set by orders of the parameters but also by the orders of the concentrations variables Xj. 
We therefore use a renormalization : 

x t = e a *Xi. (8) 

where a« are unknown powers chosen such that %{ are bounded uniformly in e (we will see later when this 
choice is possible). 

We seek for renormalization exponents such that only a few terms dominate all the others, for each 
z-th equation ^ as e — > 0. Let us denote the number of terms with minimum degree in e for i-th equation 
as irij. Naturally, 1 < to,; < Mj. After renormalization, we remove all small terms that have smaller orders 
in e as e — > 0. We can call this procedure tropical removing. The system obtained can be named tropically 
truncated system. 

Let us denote a] 3 the I th coefficient of the multi-index otij. If all to^ = 1 then we have the following 
truncated system 

^ = e^(x), F^) = P im ^ H \ (9) 
where is the index of the unique term with minimum degree in e, 



L 

1=1 

and 



Mi = Jij(i) + X! a \ a l ~ a *' ( 10 ) 



Mi < lij + X! a ? ai ~ ai jf ° r al1 3 ^ •?'(*)• ( n ) 
1=1 

If all m% = 2, in order to find possible renormalization exponents a^, it is necessary to resolve a family of 

linear programming problem. Each problem is defined by a set of pairs (j(i),k(i)) such that j(i) ^ k(i). 
We define fa by 

n n 
\ ~* , \ ik(i) 

Mi = 7ij(i) + 2^ °' ~~ a « = 7ife W + a ' _ ° 4 ( 12 ^ 
Z=l /=1 

and obtain the system of the following inequalities 

n 

Mi < lij + X! a% l ai ~ ai 0,11 3 ^ ^W- ( 13 ) 

In order to define more precisely the separation between various terms, we use the permanency concept, 
borrowed from ecology (the Lotka-Volterra model, see for instance Tak96 ). 
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Definition 3.1. The system ^ is permanent, if there are two constants C_ > and C+ > 0, a set of 
renormalization exponents asj, and a function To, such that after the renormalized variables Q satisfy 

C_ < Xi(t) < C+, for all t > Tq(x(0)) and for every i. (14) 

We assume that C± and To are uniform in (do not depend on) e as e — » 0. 

For permanent systems, we can obtain some results justifying the two procedures of tropicalization. 

Proposition 3.2. Assume that system ([6| is permanent. Let x, x be the solutions to the Cauchy problem 
for ^ and Q (or respectively, with the same initial data: 

x(0) =x(0). 

Then the difference y(t) = x(t) — x(t) satisfies the estimate 

\y(t)\ < de* exp(bt), 7 >0, (15) 

where the positive constants C\,b are uniform in e. If the original system ^ is structurally stable in the 
domain Qc-,C + — {x : C- < \x\ < C + }, then the corresponding tropical systems Q and (|5| are also 
permanent and there is an orbital topological equivalence x — h e (x) between the trajectories x{T) and x{t) 
of the corresponding Cauchy problems. The homeomorphism h e is close to the identity as e — > 0. 

The proof of the estimate (15) follows immediately by the Gronwall lemma. The second assertion 
follows directly from the definition of structural stability which means that orbits of the dynamical system 
are smoothly deformed under small perturbations. 

Permanency property is not easy to check. The following straightforward lemma gives a necessary 
condition of permanency of the system 

Lemma 3.3. Assume a tropically truncated system is permanent. Then, for each i G {1, . . . ,n}, the i-th 
equation of this system contains at least two terms. The terms should have different signs for coefficients 
Pij, i.e., one term should be a production one, while another term should be a degradation term. 

Proof. Let us suppose that to, = 1 for some i, or m.j > 1, but all terms have the same sign s. Let us 
consider this equation. Then one has, for s = 1, 

^>e»5 i (C-,C+)>0. 
dt 

Therefore, Xi(t) > 5t + 2^(0) and the system cannot be permanent. If s = —1, then 

d ^- < -e^S i (C-,C+)<0. 
Again it is clear that the system cannot be permanent. □ 



We call "tropical equilibration" , the condition in Lemma |3.3| This condition means that permanency 
is acquired only if at least two terms of different signs have the maximal order, for each equation of the 
system (J6J. 

The tropical equilibration condition can be used to determine the renormalization exponents, by the 
following algorithm. 

Step 1. For each i let us choose a pair (j(i), k(i)) such that j, k G {1, . . . , Mi} and j < k. The sign of 
the corresponding terms should be different. 

Step 2. We resolve the linear system of algebraic equations 



a l a l+ Z^ a l a '' ( 16 ) 

2 = 1 1=1 
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for ai, together with the inequalities (13). 

Notice that although that Step 2 has polynomial complexity, the tropical equilibration problem has a 
number of choices that is exponential in the number of variables at Step 1. 

Assume that, as a result of this procedure, we obtain the two terms toric system 



dxi 
~dt 



P ij± 5c a ±. 



One can expect that, in a "generic" cas^J all \ii are mutually different, namely 

A*i < A*2 < ••• < Mn-i < Mn- 



(17) 



(18) 



We can now state a sufficient condition for permanency. Let us consider the first equation (17 1 with 
i = 1 and let us denote y = xi, z = (x 2l x n ) tr . In this notation, the first equation becomes 



dy 
dt 



f{y) = bi{z)y^ -b 2 {z)y p \ h,b 2 >0, ft£R 
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(19) 



According to ( 18 ) here z(t) is a slow function of time and thus we can suppose that bi are constants (this 
step will be rendered rigorous at the end of this section, by using the concept of invariant manifold and 
methods from [Hen ). The permanency property can be then checked in an elementary way. 

Lemma 3.4. Equation \19\ has the permanence property if and only if 

Pi < h- 

For fixed z in these cases we have 

y{t, z) — > yo(z) as < oo. 



Proof. Consider the function f(y) = b\y^ x — b 2 y@ 2 . Under the condition /3i < fa, f is negative for 
sufficiently large y > 0, and positive for sufficiently small y > 0. Moreover, / has a single positive root yi 
on (0, +oo). Therefore, all the trajectories of dy/dt — f{y) tends to y\ as t — > oo and, for any S > 0, the 
interval (yi — d, y\ + S) is a trapping domain. This proves the permanency. □ 

Let us note that tropical equilibrations with permanency imply the existence of invariant manifolds. 
This allows to reduce the number of variables of the model while preserving good accuracy in the descrip- 
tion of the dynamics. The following Lemma is useful in this aspect. 

Lemma 3.5. Consider the system 

d ^ = f(y,z)=b 1 (z)i/i -6 2 (z)jA b 1 ,b 2 >6 1 >0, fteR. (20) 
dt 

**=\F(y,z), (21) 

where z £ R m , A > is a parameter and the function F enjoys the following properties. This function 
lies in an Holder class 

F e C* 1+r , r > 0, 

and the corresponding norms are uniformly bounded in Q — (0, +oo) x W, for some open domain W C R™ 1 : 

l^ 1 |c 1 +'-(fi) < C 2 - 

Supposing that multi-indices ccjj are chosen uniformly, by generic we understand almost always except for cases of 
vanishing probability, see also IGR08| 
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Assume that the condition of Lemma \3.4\ holds. We also suppose that bi are smooth functions of z for all 
z such that \z\ > 5q > 0. Assume that z €W implies \z\ > So- 
Let yx{z) be the unique solution of f{y,z) = 0. 



Then, for sufficiently small A < Aq(C2, &i, b%, /?2> Sq, ^1) equations \20) , (22) have a locally invariant 
and locally attracting manifold 

y = Y(z,X), YeC 1+r (W), (22) 

and Y has the asymptotics 

Y(z,\)= yi (z) + Y, YeC 1+r (W), (23) 

where 

\Y(z,X)\ c i +r(w) <C S \ S , s>0. (24) 



Proof. This lemma can be derived from Theorem 9.1.1 from ([Hen], Ch. 9). □ 



The generic situation described by the conditions ( 18 ) leads to trivial "chain-like" relaxation towards 
a point attractor, provided that we have permanency at each step. More precisely, all the variables have 
separate timescales and dissipative dynamics. The fastest variable relaxes first, then the second fastest 
one, and so forth, the chain of relaxations leading to a steady state. 

The following theorem describes a less trivial situation, when some timescales are not totally separated 
and limit cycles are possible. 

Theorem 3.6. Assume /ii < /X2 < ■•■ < ^ Mn holds. 

i) If the procedure, described above, leads to the permanency property at each step, where i — 1, 2, n— 



2, and if the successive application of the lemma 3.5 for the tropically truncated toric system jjlly uniquely 
ines the locally invariant smooth manifold 

Xi = Xi{x n -i,x n ), XiEC 1+r , i = l,...,n-2, r > 0, (25) 



as the unique stable hyperbolic equilibrium of the tropically truncated system (17 1. Then, the original 



system has an invariant manifold close to ( 25 ) 



Xi = Xi(x n _x,x n ) + <f>(x n _ 1 ,x n ,e), i = l,...,n-2. (26) 
where the corrections 4>i satisfy 

|&( V) e)|c*+r->0 (e->0). 

ii) If the procedure, described above, leads to the permanency property at each step, where z = 1,2,..., n— 
2, and the last two equations of the tropically truncated system have a globally attracting hyperbolic rest 
point or globally attracting hyperbolic limit cycle, then the tropically truncated system is permanent and 
has an attractor of the same type. Moreover, for sufficiently small e the initial system also is permanent 
for initial data from some appropriate domain W^^^a md has an analogous attracting hyperbolic rest point 
(limit cycle) close to the attractor of the truncated system. 

Hi) If the rest point (cycle) is not globally attracting, then we can say nothing on permanency but, for 
sufficiently small e, the initial system still has an analogous attracting hyperbolic rest point (limit cycle) 
close to the attractor of truncated system and the same topological structure. 



Proof, i) This follows from Lemma 3.5 which can be applied inductively, step by step. 

ii) Suppose that the tropically truncated system (TTS) has a globally attracting compact invariant set 
A. Let TI be an open neighborhood of this set. We can choose this neighborhood as a box that contains 
A. Then, for all initial data x(0), the corresponding trajectory x(t),x(0) lies in TI for all t > Tq(xo,TI). 
Therefore, our TTS is permanent. Here we do not use the fact that the cycle (rest point) is hyperbolic. 
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Permanency of the initial system follows from hyperbolicity of A. Hyperbolic sets are persistent 
(structurally stable Ruc89j). Since this set is globally attracting, all TTS is structurally stable (as a 
dynamical system). This implies that the initial system has a hyperbolic attractor close to A, since initial 
system is a small perturbation of TTC in 77. 

iii) If the set A is only locally attracting, the last assertion of the Theorem follows from persistency 
of hyperbolic sets. □ 



Remark. Theorem |3.6| implicitly supposes that all fast variables Xi, i = 1, 2, n — 2 can be expressed 
as functions of the remaining slow variables x n —i,x n . It does not consider the situation when the suc- 



cessive application of the lemma 3.5 leads to degenerate equilibria. This situation typically occurs when 
the tropically truncated system has conservation laws, i.e. linear combinations of the fast variables are 
invariant with respect to the truncated fast dynamics. This case, asking for variable aggregation and new 
slow variables will be discussed in detail elsewhere. 

4 Geometry of tropical equilibrations 

In this section we provide a geometrical interpretation of tropical equilibrations. We consider networks of 
biochemical reactions with mass action kinetic laws. Each reaction between reagents A t is defined as 

i k 

The stoichiometric vectors ctj € N™, f3j £ N n have coordinates ctji and (3jk and define which species 
are consumed and produced by the reaction j and in which quantities. The mass action law means that 
reaction rates are monomial functions and read 

R j (x) = k j x a =. (27) 

where kj > are kinetic constants. The network dynamics is described as follows 

dx 



dt 



(28) 



After parameters and variables rescaling, kj = kje 7j , x = xe a we obtain 

dxi 
lit 

where 



(£ t ftok j (p ji -a ji )ib a *)e- a *> (29) 



fij = 7j+ < a, a.j > . (30) 



Definition 4.1. Two reactions j, j' are equilibrated on the species i iff: 

i) fij = fi f , 

ii) ((3j - aj)i(flf - otj>)i < 0, 

iii) jj,k > Hj for any reaction k ^ j,j', such that (fik — ctfe)j 7^ 0. 



Remarks. Definition 4.1 ensures the conditions of Lemma |3.3| and is thus equivalent to tropical equili- 
bration of the species i. 



According to (30 1 and Definition 4.1 the equilibrations correspond to vectors a € R n where the 
minimum in the definition of the piecewise-affine function /j(a) = minj(7,-+ < a,aj >) is attained at 
least twice. 
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Let us consider the equality fij = fj,ji. This represents the equation of a n — 1 dimensional hyperplane 
of R", orthogonal to the vector otj — ajr, 

7j+ < a, olj >= < a, ay > (31) 

For each species i, we consider the set of reactions IZi that act on this species, namely the reaction k is 
in IZi iff (/3fc — a.k)i ^ 0. The finite set IZi can be characterized by the corresponding set of stoichiometric 
vectors ctfc. 

The set of points of E™ where at least two reactions equilibrate on the species i corresponds to the 
places where the function /j is not locally affine (the minimum in the definition of fi is attained at least 
twice). 

For each species, we also define the Newton polytope Mi, that is the convex hull of the vectors k G 



TZi. The hyperplanes defined by (31 1 and corresponding to equilibrations of two reactions on the same 
species i are orthogonal to edges of the Newton polytope Mi- Mi is also the Newton polytope of the 
polynomial -Pj(x) = • kj(f3ji — aji)x a i that defines the rhs of the ordinary differential equation satisfied 
by the species i. 

We can now state the following 

Proposition 4.2. There is a bijection between the locus 71 of vectors a where the min-plus polynomial 
fi(a) is not linear and the tropical manifold of the polynomial Pi(x) that defines the rhs of the ordinary 
differential equation satisfied by the species i. The reaction equilibrations correspond to vectors a included 



in 7j but satisfying also the condition ii) of Definition 4- 

Remarks. This property can be used to put into correspondence reaction equilibrations and slow 
invariant manifolds. Indeed, if a reaction equilibration exists, this leads to a slow manifold that is close 



to some parts of the tropical manifold of Pi{x). For instance, a reaction equilibration described by (31) 
will correspond to an invariant manifold close to a hyperplane orthogonal to ay — oty. The condition 



ii) of Definition 4.1 is needed for equilibrium (the equilibrated reactions have to have opposite effects on 
the species i, one has to produce and the other has to consume the species). Without this condition, 
the dynamics would simply cross the tropical manifold with no deviation. However, the condition ii) is 
not sufficient for stability of the equilibration (permanence). A sufficient stability condition is given by 
Lemma 3.4 and reads (ay — a'j)i(f3j — otj)i > 0. 



5 Tropical approach to the permanency problem 

We have shown in the previous sections that tropical ideas can be used to simplify complex systems, 
by tropical removing. During this procedure, permanency has to be checked at intermediate steps on 



tropically truncated systems. Lemma 3.4 allows to check permanency for toric systems with separated 
time scales. We provide here another approach to permanency, that can be applied to more general 
situations. We consider only upper estimates. The lower estimates can be found in a similar way. 

Like in the preceding sections the truncated system is obtained by removing from the non-tropicalized 
system ([T]) all the terms excepting the maximum order terms. We denote the corresponding vector field 
by F tr and the truncated differential equations read: 

1 - F *'<*>- < 32 » 

Let us first formulate a Lemma. 

Lemma 5.1. Assume that non-tropicalized system ([T]) has a smooth Lyapunov function V(x) defined on 
the cone R™ such that 

dV(x(t))/dt < (33) 
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on trajectories x(t) = (x\(t), ...x n (t)) of and 

V(x) — > oo as \x\ -» oo. (34) 
Then, if x(t) is a trajectory of ([I]) such that |x(0)| < 5', then there is a constant Cq(S') such that 

\x{t)\ < Co, t> 0. (35) 



Proof. Indeed, if \x(t)\ are unbounded as t — > +oo, one has sup t>0 V(x(t)) = +oo, but (33) entails 
V(x(t)) < V(x(0)). □ 

Let us consider the tropical version ([5]). Assume that the truncated version has a strong Lyapunov 
function V tr {x). For a truncated vector field F tr this function satisfies 

< VV tr {x{t)),F tr (x{t)) >< - K \VV tr (x(t))\\F tr (x(t))l k>0 (36) 

on trajectories x(t) = (xi(t), ...x n (t)) of ([lj and 

V tr {x) -> oo as |£c| -> oo. (37) 

Here x g R" . 

Such a function can be found for some tropical versions of two component systems. For example, if 

dx/dt = k lX a y h , (38) 

dy/dt = -k 2 x a y b , (39) 
where a, b > and k\, k 2 > 0, we can define V tr by 

F ir = x + p y , (40) 

where [3k 2 > k x . Then W = (l,/3), and Q, d37} hold. 



Lemma 5.2. Assume the tropicalized system ([5]) has a smooth Lyapunov function y' r (a;) defined on the 
cone R™ such that (36 1, (37 1 hold. Assume that Fi(x) and Ff r (x) are multivariate polynomials of x such 
that deg{Fl r ) > deg{F l ), where F, = F t - Fl r . 

Then, if x(0) < 8', then there is a constant Cq(5') such that solutions of non-tropicalized system Q 
satisfy 

\x(t)\<C , t>0. (41) 

Proof. Let us compute the derivative dV tr jdt on trajectories of the initial (non-tropicalized) system. We 
have the relation 

dV tr (x(t))/dt =< VV tr {x{t)).F tr {x{t)) > + < VV tr (x(t)),F(x{t)) > . (42) 



Using the definition of strong Lyapunov functions, from (42 1 one has 

dV tr (x(t))/dt < \VV tr (x(t))\^ K \F tr (x(t))\ + \F(x(t))\). (43) 



But for large \x\ one has |F(a;)| < K,\F tr (x(t))\, because \F(x)\ = o(|F tr (a;)|), \x\ -> 00. Therefore, (|43j) 
gives then 

dV tr (x(t))/dt < 0. (44) 

This shows that |a;(t)| cannot increase to +00, and finishes the proof. □ 
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6 Application to chemical reactions kinetics 

As an application, we discuss the Michaelis-Menten mechanism of catalysed reaction. This model can be 
schematically described as : 

S + E ^ ESH P + E, 

k-i 

where S, E, ES, P represent the substrate, the enzyme, the enzyme-substrate complex and the product, 
respectively. 

The rate functions obey mass-action laws. We denote by x — [S] and y = [SE], the concentration 
of substrate and enzyme-substrate complex, respectively. The reaction mechanism has two conserved 
quantities eo = [E] + [ES], so = [S] + [ES] + [P]. Using the conservation laws we obtain the following 
reduced system: 

x' = -fcix(e - y) + k-xy, 

y' = kxx(e — y)- (ft-i + k 2 )y. (45) 



Let us consider that the initial data satisfies < y(0) < eo, < x(0) + y(0) < sq. Then, from (45) it 
follows that 

< y < eo, < x + y < s , < x. (46) 

This type of constraints are typical for reduced systems resulting from ODE systems with conservation. 
We solve now the tropical equilibration problem. Using rescaled variables x = xe ai , y = ye a2 , k\ = 
fcie 71 , k-i — fc-ie 7 " 1 , eo = eoe 7c , sq = sge 73 , (45) becomes: 

x' = -fcie e 7l+7e S + k 1 £~ <1+a2 xy + fc_ie 7 - 1+a2 " ai y, 

y' = fc 1 e e 7l+7e+ai_a2 x - Aiie 7l + ai sy - (k^- 1 + k 2 e l2 )y. (47) 

The two tropical equilibration equations for x and y read: 

7i + 7e = min(7i + a 2 , 7-i + a 2 - ai), (48) 
7i+7e+ai-a 2 = min(7i + ax, min(7_i, 72)). (49) 



We should add to these, the constraints ( 46 ) imposed by the dynamics: 

»2 > 7e, min(ai, a 2 ) > 7s- (50) 

We can distinguish between two situations. 

Let us first consider that 7-1 < 72- In this case (49) is equivalent to (48) (it can be derived from 
the latter by adding ax — ol 2 to both sides). This situation corresponds to much larger than k 2 and 
means that the enzyme-substrate complex is recycled to a much larger extent than it is transformed into 
the reaction product. We can find two solutions for the tropical equilibration problem and two different 
tropically truncated systems (TTS). 

The first solution demands large concentrations of substrate and corresponds to saturation of the 
enzyme (saturation regime): 

a x < 7-1 - 7i) fl 2 = 7e, 
x 1 = e 7l+7 = (— kie x + kxxy), 
g = e 7i+«i (kxeox - kxxy). (51) 

The second solution works for small concentrations of substrate (linear regime): 

ai > 7-1 - 7l. «2 = ai + 7e + 7i - 7-1: 
x 1 = e 7l+7e (— kxeox + k-xy), 
y> = et-i{kxe x-k-iy)- (52) 
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In order to further characterize these two functioning regimes, we consider the third variable z — (x + 
y)e~ 7s . The choice of this variable is dictated by the TTS. In general, conserved quantities of the TTS 
(total mass of fast cycles) can be slow variables of the full system [RGZN12]. If this variable is slower 
than both x and y, the regime is called quasi-equilibrium IGRZ10|lGorll| and consists in rapid exchanges 
between substrate and enzyme and a much slower transformation of the total mass [S] + [SE] into [P]. 
In both cases the equation for z reads: 

z' = -e^+^-^hy- (53) 

A sufficient condition for quasi-equilibrium (ensuring both 72 + a 2 — 7 S > max (71 + 7 e ,7i + »i) and 
72 + ct2 — 7s > max (71 + j e , 7-1) in the first and second of the cases above, respectively) is 72 > 71 + j s - 
The second situation is when 7_i > 72. This case leads to negligible recycling of the enzyme-substrate 
complex that is rapidly transformed into reaction product. Quasi-equilibrium is no longer possible, but 
we have another interesting equilibration corresponding to fast consumption of one of the variables. The 
QSS variable is necessarily equilibrated and fast. The remaining variable is slow. This corresponds to the 
well known quasi-steady state (QSS) regime of the Michaelis-Menten mechanism, first discussed by Briggs 
and Haldane |GRZ10( IGorlR IRGZN12] . 



In this case (48 1,(49) are no longer equivalent: 



7i + 7e = min(7i + a 2 ,7-i + a 2 - ai), (54) 
7i+7e+ai-a 2 = min(7i + ax, 72). (55) 

We obtain four solutions to the tropical equilibration problem and four different truncated systems. In 
three of these solutions, only one variable is equilibrated (see Table 1). The solutions 1 and 4 correspond 
to rapid complex consumption in saturated and linear regimes, respectively. It is the case discussed by 
Briggs and Haldane. The solutions 2 and 3 correspond to very small concentrations of the substrate. 



Table 1: Tropical equilibrations of the Michaelis-Menten model with negligible recycling 72 < 7-1. All these 
equilibrations have a geometrical interpretation illustrated in Figure 1. 



No 


Condition 


Truncated system 


Regime 


1 


a\ < 72 - 7i 
a 2 = 7e 


x 1 = e 7l+7<! (— fcieox + k\xy) 
y* = e 7i+ai (kieox — kixy) 


y QSS if 

a\ < 7 e 


2 


72 - 71 < ai < 7_i - 71 
a-2 = 7e 


x 1 = eTi+T<= (—ki&ox + k\xy) 
y' = -e l2 k 2 y 


x QSS if 

72 > 71 + 7e 


3 


ai > 7_i - 71 

a 2 = ai + 7 e + 71 - 7_ x 


x 1 = e 7l+7e (— fcieox + k-iy) 
y> = - e 72^ 2 y 


x QSS if 

72 > 71 + 7e 


4 


ai > 72 - 7i 

a-2 = ai + 7 e + 71 - 72 


x 1 = — e 7l+7o fcieoa; 
y' = e j2 (kieox - k 2 y) 


y QSS if 

72 < 71 + 7e 



7 Conclusion 

Tropical analysis provides useful tools for understanding the dynamics of biochemical networks. In this 
paper we have studied the simple example of an enzymatic reaction, but some other applications have 
been discussed elsewhere, see [NGVR121 |SCF + Q9] . We have shown that depending on the values of the 
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Figure 1: Newton polygon and tropical manifolds for x (in blue) and for y (in red) variables of the Michaelis-Menten 
model. If 7_i < 72 the two manifolds coincide in the limit e — > 0, meaning that both variables are equilibrated. 
The vertices of the Newton polygons correspond to monomial terms in the ODEs (different vertex shapes mean 
different signs of the monomials). Only two edges of the Newton polygon relates vertices of opposite signs, which 
means that there are two equilibrations; these correspond to the branches OA3 = O'As and OA\ = 0'A[ of the 
tropical manifolds. If 7-1 > 72, the two tropical manifolds for x and for y share a common half-line, but no 
longer coincide. This leads to four possible equilibrations as in Table 1: O'A'l (solution 1), OO' (solution 2), OA3 
(solution 3), and O' A' 3 (solution 4). 



parameters and concentrations, biochemical networks with multiple time scales can have several asymptotic 
regimes. During such regimes, the dynamics can be approximated by truncated systems obtained by 
tropicalization of the ordinary differential equations describing the chemical kinetics. Tropical geometry 
can guide the construction of such truncated systems. An important step in this construction is the 
calculation of tropical equilibrations leading to slow invariant manifolds. We showed that there is one to 
one correspondence between tropical equilibrations and well defined parts of the tropical manifolds of the 
polynomials defining the ordinary differential equations. In the future, effective algorithms will be needed 
for the tropical equilibration problem, whose complexity is exponential in the number of variables. This 
will be essential for large scale applications in systems biology. Also, our methods will be generalized to 
include the case when the tropically truncated fast subsystem has conservation laws, when aggregated 
slow variables are needed. 
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